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.^ I A theory of the quasidilaton is an extension of massive gravity by a scalar field, nonlinear ly realizing a 

fvq \ certain new global symmetry of the Lagrangian. It has been shown that unlike pure massive gravity, this 
theory does admit homogeneous and isotropic spatially flat solutions. Among the latter, selfaccelerated 
solutions attract a special attention. Previous studies of perturbations, performed in the decoupling limit, 
f^ revealed one healthy scalar mode, while the second relevant scalar was not captured in that limit. Here we 

^ ' study full cosmological perturbations above the simplest selfaccelerated background. We show that the 
fluctuations of a mixed state of the quasidilaton and the helicity-0 graviton necessarily have a negative 
kinetic term at short distances, making this background unphysical. In addition, these cosmologies exhibit 
an order one sensitivity to higher dimensional terms suppressed by an energy scale that is parametrically 
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CSJ ' higher than the strong coupling scale of the quasidilaton effective theory: such terms include Galileons, 



Goldstone-like selflnteractions and derivatives of the quasidilaton coupled to curvature, none of which 



.^ ' introduce extra Ostrogradsky states. As one consequence, cosmology at the Hubble distances for this 
^^ ' particular class of solutions depends on an unknown extension of the quasidilaton below its strong coupling 
distance scale. We note that non-FRW solutions that are similar to those of pure massive gravity should 



not necessarily suffer from these pathologies. 



I. INTRODUCTION AND SUMMARY 

The theoretical robustness and remarkable accuracy of General Relativity (GR) in describing 
gravitational interactions in a wide range of distance scales makes it one of the most successful 
physical theories of all times. Nevertheless, the late-time cosmic acceleration, and the cosmological 
constant problems might be pointing towards a certain missing ingredient in the GR picture of 
gravity, calling for its modifications at distances of the order of the present-day Hubble scale. Mas- 
sive gravity is perhaps the most conservative of such modifications. Whether or not an interacting 
graviton can have a nonzero mass without violating the theoretical consistency of GR has been a 
long standing problem originating from the work of Fierz and Pauli [1], who proposed a unique 
linear theory, propagating the five degrees of freedom of a massive spin-2 state. Generalizing the 
theory to the nonlinear level, however, has proven to be problematic, due to the emergence of an 
extra ghost-like state [2]; this 6th degree of freedom is referred to as the Boulware-Deser (BD) 
ghost. 

Recently, an explicit construction has been proposed in Ref. [3], that gave an order-by-order 
Lagrangian for a single massive spin-2 state free of the BD ghost in a particular limit. In Ref. [4] 
the Lagrangian of [3] was resummed into a diffeomorphism invariant nonlinear theory, which was 
proposed as a ghost-free candidate for massive GR. The presence of the hamiltonian constraint, 
necessary for projecting out the BD ghost was shown perturbatively, up to the quartic order in 
Ref. [4]. The full proof of the ghost-freedom to all orders was given in [5]. Subsequently, the 
theory has also been shown to be BD ghost-free in the Lagrangian [6], as well as in the vierbein 
formalisms [7, 8]. 

An important property of the ghost-free theories of massive gravity is the existence of self- 
accelerated solutions with the Hubble scale of order of the graviton mass, as pointed out in Ref. 
[9]. Since then, these theories have been widely studied in the context of cosmology [10-17], spher- 
ically symmetric solutions and black holes [18-23], as well as their quantum consistency [24, 25]. 
One interesting feature of cosmology in these theories is the absence of strictly homogeneous and 
isotropic spatially flat and closed backgrounds, while the obtained inhomogeneous solutions may 
still well approximate the observed world [11]. In particular, there are self-accelerated solutions 
for which the metric can be brought to the standard flat FRW form at the expense of having 
inhomogeneities in the Stiickelberg fields [11, 14, 26]. The spatially open FRW solutions have been 
found in [12]. However, the above types of solutions generically suffer from nonlinear instabilities^ 
[13, 27, 28]. 

Recently, it has been shown that an extension of massive GR by a scalar, the quasidilaton, 
nonlinearly realizing a certain new global symmetry, naturally reintroduces homogeneous and 
isotropic, spatially fiat solutions, for certain values of the parameter characterizing its coupling to 
gravity [29]. Such solutions include selfaccelerated backgrounds that seem to be consistent with 
all immediate cosmological/astrophysical tests [29]. An important property, that automatically 



^ While the theories at hand are free of tlie BD ghost, one or more degrees of freedom out of the physical 5, can 
flip the sign of their kinetic terms on a nontrivial background, rendering it unstable. This is what we refer to as 
"ghost instability" here. 



follows from the construction of the theory, is its technical naturalness. Indeed, the form of the 
interactions of the quasidilaton is protected by the new global symmetry of the theory. Moreover, 
in the decoupling limit, quasidilaton massive gravity (QMG) acquires yet another, enhanced global 
symmetry, reducing to a theory of two Galileons [30] interacting with a tensor field [3]; all these 
terms exhibit remarkable non-renormalization properties [24, 31, 32]. While the theory certainly 
stands out in this respect, it is crucial to address the question of stability of the spectrum of 
perturbations on such selfaccelerated backgrounds. 

The study of perturbations on the simplest selfaccelerated de Sitter (dS) solutions of QMG has 
been initiated in [29], in a certain high-energy (decoupling) limit of the theory. It has been shown 
however, that the decoupling limit analysis only captures one of the two scalar degrees of freedom; 
whether or not the second scalar propagates and is stable, remained an open question. In this 
Letter we perform the full analysis, and show that the second scalar is dynamical. However, we 
find that the kinetic term of this mode necessarily has a negative sign at short distances, rendering 
these simplest dS backgrounds unstable. 

We note, in addition, that these backgrounds exhibit sensitivity to an unknown short distance 
physics: the solution is fully modified by adding the Galileon and/or Goldstone-like selfinteractions, 
as well as certain derivative couplings to curvature tensors, which are consistent with all the 
symmetries of the Lagrangian, and do not introduce any Ostrogradsky ghost states. Some of 
these higher dimensional terms are suppressed by a scale that is parametrically higher that the 
strong coupling scale of QMG. The order one sensitivity to such terms makes the viability of these 
particular backgrounds questionable, even though the selfaccelerated solutions do generically exist 
even when these terms are included. In particular, even in the full quantum theory, where such 
terms are expected to be generated, one should anticipate selfaccelerated backgrounds with the 
Hubble parameter H of the order of the graviton mass; however, as pointed out above, the precise 
nature and detailed properties of the spectrum of perturbations on these solutions are fully sensitive 
to an unknown extension of the quasidilaton below its strong coupling distance scale. 

While such a UV-sensitivity is a characteristic feature of the particular class of flat FRW 
solutions considered below and, as we will see, stems from the large rate of variation of the 
quasidilaton field on these solutions, there should exist other solutions that closely resemble the 
inhomogeneous and/or anisotropic cosmologies of pure massive GR. In particular, for sufficiently 
large values of the parameter w, that controls the strength of the coupling of the quasidilaton to 
the rest of the fields, one expects the existence of inhomogeneous solutions that can recover the 
standard early cosmology to a great accuracy [11] . 

Moreover, even for the values of u: discussed in this work, < cj < 6, one may explore the 
existence of UV-insensitive inhomogeneous solutions, characterized by a small (or zero) expectation 
value of the quasidilaton, relying on the Vainshtein mechanism [33] in the cosmological setup [11] 
for the other fields. This, however, is outside of the scope of the present work. 

II. PERTURBATIONS ON THE SELF-ACCELERATED BACKGROUND 

The framework: The theory we wish to consider is based on the recently formulated class of 
models of massive gravity, free of the Boulware-Deser ghost [3, 4]. We introduce a special scalar a 



that gives rise to a certain new global symmetry of the Lagrangian. The symmetry transformation 
involves the scalar itself, and the Stiickelberg fields 0", a = 0, 1, 2, 3, that are necessary if one 
wishes to work with a diffeomorphism-invariant action for massive GR. These four fields are scalars 
w.r.t. diffeomorphisms, but do transform under the Poincare group of the internal space of ^'^'s, 
as emphasized in [25, 34]. The new global symmetry that we use as a building principle for the 
gravitational action involving the scalar a^ is realized as follows: 

a ^ CT - a Mpi , 0" ^ e"0" , (1) 

where a is an arbitrary symmetry transformation parameter. The rest of the fields in the Einstein 
frame^, and the physical coordinates x^, do not transform. This symmetry fixes uniquely, modulo 
derivative interactions, an extension of massive GR by the a field; in particular, one consequence 
of (1) is minimal coupling of matter to gravity in the Einstein frame (unlike Brans-Dicke theories 
for instance, which have matter coupled minimally to gravity in the Jordan frame). 

The most general Lagrangian invariant under (1), and excluding possible ghost-free derivative 
interactions^ of a, is written in the Einstein frame as follows [29] 



+ ?,Ml,m'' j d^x P e^'^/^^P' ^-det d^rdu'Pa + j d'x y/^C.^{g,^,ij) . 



(2) 



Here we have defined K,^^ = ^u ^ e^/^^pi ^ygt^°'^a(|)°'^u(f)^f]ab , where r]ab is a matrix, that numerically 
coincides with the flat metric (we will use the mostly plus signature), and Ui are its specific 
antisymmetric polynomials 

^2 = 2£^„. e'^^-Xt'^1^ (3a) 

U, = e^^.e^^'-t^X'^^Vs (3b) 

U, = e^^,pe''^'''iC^^X''pt\iCP, . (3c) 

The dimensionless constants 03^4, /3, cu, represent the four free parameters of the theory (the mass 
is another free parameter, but for phenomenological reasons, we set it to be of order of the present 
Hubble scale.) We will refer to the first two terms of (2) as the 'gravity sector', while the last 
one specifies the coupling of matter fields ip to it, which, in order to respect the symmetry (1), 
is given by the regular general- relativistic minimal coupling to the Einstein frame metric g^^. As 
emphasized above, the defining property of the gravity sector is its invariance under the global 
one-parameter group (1), which is related to (nonlinearly realized) scale transformations and is 
explicitly broken by couplings to matter. 



^ We define the Einstein frame in the standard way - the one for which the kinetic term for the graviton has the 
usual Einstein-Hilber form; Jordan frame on the other hand wiU feature a kinetic mixing between the scalar a 
and the graviton. 

^ We will return to the derivative interactions in the last section. 



Background cosmology: The theory at hand has been extensively studied in Ref. [29]. In 
particular, it has been shown to be BD ghost-free at the full non-linear level; moreover, it has 
been found that in the decoupling limit it reduces to a bi-galileon theory, thus inheriting all of 
the properties, characteristic of the galileon field theories [30]. This includes ghost-freedom, a 
successful implementation of the Vainshtein screening, non-renormalization, etc. (see also [35, 36] 
for models of the early universe, based on galileons.) Most importantly for the present work 
however, quasidilaton extended massive gravity has been shown to admit strictly homogeneous 
and isotropic flat FRW cosmologies - in sharp contrast to the pure massive gravity, as well as to 
its previously studied extensions. The absence of the standard cosmologies in the latter theories is 
directly linked to ghost-freedom. Indeed, the same constraint that removes the BD ghost from the 
spectrum, severely constrains the homogeneous and isotropic evolution of the universe, essentially 
restricting the scale factor of such a universe to be static, d = 0. Standard FRW evolution 
however becomes possible in QMG, due to the presence of an extra scalar field, whose cosmological 
variation allows to lift the above restriction. We briefly describe how this works in the appendix, 
concentrating on the homogeneous and isotropic field configurations in the Einstein frame theory 
and considering the asymptotic dynamics, when the matter/radiation has redshifted away^. In 
particular, in terms of the three free parameters of the theory (we assume the graviton mass is of 
order of the present Hubble scale), the Hubble constant is given by 



H' 



"^^ ( i("3 + 4a4)c^ - (1 + |«3 + 3a4)c^ + (3 + fag + 3a4)c - (2 + ag + a^) j 



(4) 
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where c denotes a specific combination of the three parameters a^^^ and /3, which solves the quartic 
equation, given in (A2). For the particular case of /3 = 0, the two interesting solutions are given 
by 

3a3 + 8^4 ± \/9al — GAa^ , . 

The de Sitter solutions exist in the latter case as far as the following conditions are met [29] 

q;2 

ag ^ 0, < 04 < — , < a; < 6 . (6) 

8 

The last condition, < w < 6, is necessary for the well-definiteness of the solution at hand, 
regardless of whether (3 is zero or not. 

Perturbations: The study of perturbations on the above self-accelerated background has been 
initiated in [29] , where it has been shown that the conditions for the existence of the dS solutions 
also guarantee that ghost instabilities are absent for both the tensor and vector perturbations (see 
also [37] for a discussion of stability in a general Lorentz-breaking massive gravity on cosmological 



^ We refer the reader to [291 for more details. 



backgrounds.) The remaining open question concerns the ghost-freedom of the scalar sector. 

As we will see shortly, the situation is complicated by an intricate nature of the scalar pertur- 
bations. There can be up to two propagating scalar degrees of freedom on a generic background 
(the helicity-zero graviton and the quasidilaton) . In massive gravity and its closely related theo- 
ries however, cosmological backgrounds commonly happen to be degenerate, with some degrees of 
freedom losing their dynamics, or becoming infinitely strongly coupled. In QMG, one can show 
that the decoupling limit treatment only captures one mode [29], without being able to say any- 
thing about the other; one therefore has to resort to the full treatment of the scalar perturbation 
Lagrangian, which we perform in the rest of the present section. The resulting analysis is rather 
technical. We will suppress lengthy expressions whenever possible, providing instead a detailed 
account of the procedure used and results obtained. 

Our starting point is the perturbation Lagrangian given in (A9) - (A12) in the unitary-like 
gauge. Various conventions and definitions used are given in the appendix; in particular, we will 
work with the conformal time r, with a prime denoting differentiation with respect to it. The 
action, charecterized by the parameters 7^ explicitly given in (A12), involves the perturbation of 
the quasidilaton (, coupled to the metric perturbations /iqo, ^oj and hij. 

It is relatively easy to deal with the tensor and vector modes (see [29, 37].) In particular, one 
can show, that as far as 75 > 0, there are no ghost instabilities in either of the two sectors. For the 
particular case of /3 = 0, this always holds as far as the dS existence conditions (6) are satisfied. 
Gradient instabilities on the other hand never occur in the tensor sector, while their absence for 
the vector modes requires 

76 < , (7) 

which is also the condition for the absence of tachyonic instability for the tensor modes. 

Turning to the discussion of the scalar sector, we present the essential steps of the calculation, 
along with the basic results. 

With the following decomposition 

did- 
/loo = V', hoi = diV, hij = -^p + 6ij(f) , (8) 

the Lagrangian for scalar perturbations in the spatial momentum space^ can be written as follows 
(to simphfy expressions, we will make use of the following relations, •Jq = —77, 74 = —71/6, that 



^ To avoid complicating the presentation, we will not change the notation for the spatially Fourier-transformed 
fields. 



hold for the most general choice of the parameters) 



a^ 



^-T 



^(0' + aHi^f - lk20(2^ - 0) - (0' + aHi^){2\^v + p') + uj{C - k'C') 



+ ujaH{^l) + 30 + p)C' - 2LoaHk\C + a^7iV^C + a^72(30 + p)C + a^73^^ + 0^74(80 + p)^ 
+ a^-f^k^v'^ - 20^76(302 + 20p) + a'^-fsC^ 

An immediate observation is that ip and v enter without time derivatives. Integrating out these 
fields, one obtains a fairly complicated effective Lagrangian £3 for p, and (, the exact form of 
which we will not reproduce here. An important property of this Lagrangian however, is that the 
two-derivative kinetic matrix K3, defined for the fields $ = (p, 0, Q"^ by 

£3 = $'^i^3$' + . . . , (9) 

is degenerate, having one zero eigenvalue with the corresponding eigenstate being pi = 20+C. This 
combination is thus non-dynamical and will have to be integrated out in order to obtain the effective 
action for the two remaining dynamical states. Of course, the vanishing of the determinant for the 
kinetic matrix K3 should be no surprise (and is a nice consistency check of the computation), since 
the BD ghost-freedom of the theory guarantees that there are no more than two scalar degrees of 
freedom propagating on an arbitrary background. We take the latter two to be 

p2 = -^(p' + C, P3=P ■ 

Solving the pi - equation of motion and substituting the solution back into the action yields the 
final Lagrangian £2, that determines the dynamics of the two propagating scalars. In the short- 
wavelenth limit, (|k| ^ m), the kinetic matrix K2 for the fields p2 and p3 drastically simplifies: 
it is degenerate, with the non-zero eigenvalue being Z^^ = a^uj jl for the eigenstate p2, while the 
kinetic term for p3 vanishes in this limit®. In order to capture the dynamics of p3 therefore, we 
have to go to higher order in the |k|~^ expansion. The easiest way to do this is to deal with the 
determinant of the kinetic matrix, which has the following expansion 

In order to obtain the latter equation, the following relations have been used 

^1 = -7^ + 67175 + Qwi/S = 0, ^2 = 672 + 2477 + 78 - 9wif ^ = , (11) 

which can be obtained on the basis of the constraint equation (A2), and drastically simplify the 
expressions at hand. The leading contribution to the second eigenvalue of the kinetic matrix in 



This explains the absence of the second scalar mode in the decoupling limit, as found in [29]. 



the UV is therefore 



^vv _ a^(w - 6)7? 



3 48a;k4 

Unfortunately, the condition for the existence of a well-defined de Sitter solution, a; < 6, constrains 
this quantity to be negative, thus signalling a ghost instability for the short wavelength modes. 
The full expression for the determinant, using the relations (11), can also be written in a simple 
form 

det K, = , ^^(6-"h? . (12) 

One can see that in going from the UV to the IR, the determinant changes sign exactly once, being 
manifestly positive for the long-wavelength modes. 

UV dispersion relations: As we will see momentarily, the vanishing of the 1/k^ piece in 
the expansion of the determinant (proportional to ^i) is not accidental: it is required in order to 
obtain the regular dispersion relations E'^ ~ k^ for the short wavelength modes^. The kinetic part 
of the Lagrangian in the k^ m limit can be written as follows: 

where ellipses denote possible corrections of higher order in the short-distance expansion. There 
is a correction to one of the eigenstates of the kinetic matrix, 

7i(a;-6)a^ ,^ ^, 

"'="'+ 12k J "■'' ("' 

which is retained above, since it gives an 0{1) effect in the UV, as will become clear below (this 
can easily be seen by canonically normalizing the pa field.) 

As a next step, we have a look at the one-derivative part of the action. There are four possible 
types of terms: P2P2' P3P35 P2P3 and P2P3- The first two of these are total derivatives and, upon 
integration by parts, can be reduced to non-derivative mass/gradient terms (we should keep in 
mind that these operators come with coefficients that depend on the conformal time through a{T), 
thus giving nonzero contributions.) It is convenient to arrange the rest of the one-derivative terms 
into combinations P2P3 ± P^P'i- The first of these is again a total derivative and therefore only 
contributes to the non- derivative part, while the hermitian antisymmetric combination starts from 



^ Of course, the fact that the leading UV contribution to the kinetic term of the pz field goes as k "^ stems from 
the specific normalization of the metric perturbation a in (8). 



the order 1/k^ and enters into the Lagrangian as follows ^ 

^r ^ § (P2P^ - pM . (15) 

At the level of zero time derivatives, there are three types of terms: p\, p2p3 and p^. Keeping track 
of all partial integrations done up to this stage, and leaving only the leading contributions in the 
short-distance expansion, the relevant terms are given as follows 

rUV ^ 1 2 2 I I '^OS 2 I n r\ 

£2 =) Coik P2 + C02P2P3 + -^Ps + ... . (16) 

In order to extract the physical scalar spectrum on the given de Sitter background, we will have to 
work in terms of the eigenstates of the kinetic matrix p2 and ps, canonically normalized^ as follows. 



P2 — > P2I y 'Z'2^ , Pi — > pzl yZ^'^ . This procedure puts the kinetic term in the form, invariant 
under orthogonal transformations of the two dynamical fields, that one can use to diagonalize the 
mass/gradient matrix. For perturbations of frequency E, exceeding the characteristic time of the 
Hubble expansion, it is convenient to Fourier transform in time as well, ' ~ E'. It is easy to see 
now that the kinetic terms are of order E'^ at high energies, one-derivative terms are of order mE, 
while the gradient terms are ~ k^. This shows that the dispersion relations are of the regular 
£" ~ k form in the UV, so that the one-derivative contribution can be completely neglected. 

III. UV SENSITIVITY 

We have shown above, that the scalar perturbation spectrum on the homogeneous and isotropic 
de Sitter solutions obtained in quasidilaton extended massive gravity in [29], for the most general 
choice of the free parameters, features all dynamical modes present in the theory. Moreover, one 
combination of scalars is necessarily a ghost with the regular dispersion relation E'^ ~ k^ at short 
distances, rendering these backgrounds unstable. 

However, the theory considered above is not the most general one, consistent with the global 
symmetry (1). Indeed, while we have certainly included all possible non-derivative interactions of 
the quasidilaton field, there still are ghost-free derivative interactions that may be added to the 
quasidilaton Lagrangian. For instance, one could add the covariantized cubic Galileon (see [30] for 
general Galileons, and [38] for their covariantization) 

v^ A3 • (17) 

This term is invariant w.r.t. shifts of a by a constant, and hence preserves the symmetry (1). 



^ There is a seemingly accidental cancellation in that the piece of order k" only enters with the total derivative 
combination P2P3+P2P3J but not with the orthogonal one. Again, this has to be the case if the dispersion relation 
is to be quadratic in momenta in the UV. 

^ Note, that the normalization factors Z explicitly depend on time through the scale factor. Canonical normalization 
will thcreforrc produce extra one and zero-derivative terms. These will however be subleading at short distances. 



Moreover, on an arbitratry background it gives rise to the second order equation of motion. Note 
that here we chose the suppression scale of this operator to be A3 = M-piim?, since this is the strong 
couphng scale of the quasidilaton theory on the Minkowski space [29]. At this scale the theory 
needs a UV extension, and adding such a term would in general be motivated in the full quantum 
theory. 

Can this and other Galileon terms affect the conclusions of our previous section? At first sight 
the answer seems negative, since the selfaccelerated solution we are discussing is characterized by 
the curvature scales of order l/m^ that are much greater than A^^. However, this naive expectation 
is questionable: for the homogeneous and isotropic cosmology with these new terms we still have 
the relation 

a ~ MpiH, 

which is not modified by the inclusion of the Galileons. Due to this scaling, the effects of the 
Galileon terms cannot in general be ignored in our considerations. 

Furthermore, along the same lines, one could add to the quasidilaton Lagrangian the symmetry- 
preserving Goldstone-like self-interactions of the a field: 

V9 74^134 > U»J 

where the natural scale for these operators is A2 = (Mpi?Ti)^/^. Note that A2 >> A3. Therefore, 
these operators are suppressed by the scale that is higher than the flat-space strong coupling scale 
of QMG. Naively, such terms should be irrelevant at scales 1/m. Nevertheless, due to the relation, 
a ~ Mp\H ~ A2, these terms do modify the selfaccelerated background as well as properties of 
perturbations on it. To see this more explicitly, let us for simplicity concentrate on the term with 
n = 2. Since it does not include the Stiickelberg scalars, the constraint equation in (A2) is not 
modified, giving a ~ MpiH, while the Friedmann equation receives a correction proportional to 
H'^. Restoring all the scales, the Friedmann equation takes the following schematic form: 



yiH^ + y2—^ = y-m\ (19) 



where the coefficients yi are some functions of the dimensionless parameters of the theory. One 
can therefore anticipate the de Sitter solutions with H ^ m, just as in the theory without the 
derivative a self-interactions. Moreover, the structure of perturbations on this solution will also 
get modified, since the terms (18) generically contribute to all pieces in the quadratic perturbation 
Lagrangian (AlO). 

Moreover, even if the tree-level theory is taken to be the one considered above, quantum loops 
should generate the operators of the form (18). Naively, the magnitude of these operators in (18) 
seems to be smaller than what one would expect from knowing that the cutoff of the theory is A3. 
However, the fact that they should be suppressed by a scale higher than A3, stems from the special 
nature of the Minkowski space decoupling limit of the theory. Indeed, as noted in [29], in the 
decoupling limit the theory acquires an enhanced galilean symmetry, under which the quasidilaton 
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shifts as 

a ^a + b^x^ . (20) 

Goldstone-type interactions are not invariant under this transformation, and would not be gener- 
ated in the decouphng hmit, if absent from the tree-level theory. This means that whatever the 
correction is in the full theory, it should not survive taking the decoupling limit, implying that A2 
sould be taken as the suppressions scale for this operator. 

The Galileon and Goldstone-like interactions do not exhaust the full list of all possible ghost- 
free extensions of the model. Indeed, there are certain derivative couplings to the curvature 
tensors [20, 39], that represent a covariantization of the ghost-free scalar tensor interactions of the 
decoupling limit of massive OR [3]. In particular, one can extend the Lagrangian by the following 
interaction 

Mpim v^ -^ — -g , (21) 



Ai 



as well as a higher-order term 



A? 



M^,m^ ^^-^^ \y^ , (22) 



where L^avp is a fairly complicated combination of the Riemann tensor and its contractions, that 
can be found in [20, 39]. Both of these terms (as well as their arbitrary functions) are invariant 
under the quasidilatations (1) and are as important on the self- accelerated FRW backgrounds, as 
the rest of the terms considered above; they would be expected to affect the analysis at order 
one. Moreover, as already emphasized above, these terms lead to the Ostrogradsky ghost-free 
scalar-tensor interactions of a with the metric perturbation in the Minkowski space decoupling 
limit, retaining the remarkable non-renormalization properties of the theory. 

It is important to stress, that while the very existence of selfaccelerated solutions with the 
Hubble parameter of order of the graviton mass might also be the property of the full quantum 
theory, the nature of these backgrounds and their perturbations seems to crucially depend on an 
unknown UV extension of the theory at and above the A3 energy scale. Knowledge of a putative 
UV theory with the quasidilaton, valid all the way up to the Planck scale, becomes a necessity for 
the full understanding of homoegenous and isotropic flat FRW backgrounds in this theory. 

It is clear that for the fiat FRW solutions at hand, due to the large time variation of the 
quasidilaton field, a ~ A|, both the late time {H r^ m), as well as the early {H ^ m) cosmology 
are sensitive to the unknown UV dynamics of the theory. This is unlike the pure massive gravity 
theories, for which all of the previously obtained solutions are well within the effective field theory 
domain, possible due to the parametric separation of the cosmological and the UV scalaes, m <^ A3. 
Now, in the limit w — )■ 00, the a field is expected to decouple and the dynamics of the theory, along 
with its solutions, should be dominated by the pure massive OR [29]. In particular, for sufficiently 
large a;, one expects the existence in these theories of the inhomogeneous solutions of massive GR, 
that recover the standard early cosmology due to the Vainshtein mechanism [11]. 
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Appendix A 

The most general ansatz for flat, homogeneous and isotropic solutions is given as follows: 

ds"^ = -N^{t)dt^ + a\t)dx\ (tP = f{t), (f)' = x\ a = a{t). (Al) 

One can substitute this into (2) to obtain the minisuperspace action^°. Varying the resulting 
expression w.r.t. / (that plays the role of a Lagrange multiplier in the action) yields the modified 
constraint, which can be solved by 

<T/Afpi 



ca, c ( 1 + -tta + ^4 - ( 1 + -as + Sa^ j c + - (as + Aa^) c - (04 + /3)c j = (A2) 

The Friedmann equation is obtained by varying with respect to TV and subsequently setting A^ — ;■ 1. 
This yields a de Sitter metric with the Hubble scale H of order of the graviton mass m. 

"^^ ( i(«3 + 4a4)c^ - (1 + i«3 + Sa4)c^ + (3 + fag + 3a4)c - (2 + ag + 04) ) 

tt2 \ J (A3) 

^ 6 

Finally, the Lagrange multiplier / = 0° can be obtained from the a equation of motion and is 
given as follows 

iD - f dt - ^ u H"^ , , ,, 

(f) = c -—-, c = const = 1 H -, (A4) 

J a{t) n m^ 



where 

3 \ . Z'. 3 „ \ 3 



K = C 



3(1 + ^as + 04 j -2(1 + ^^3 + 3a4 | c + j (0^3 + 4^4) c^ 



(A5) 



To study the perturbations, it will be more convenient to work in terms of conformal time 
T, transforming to an "almost unitary" gauge in which the background metric is g ^ = a'^i^r)'!] ^ 



^° Note that wc have retained the lapse N{t) in the action despite the fact that by time reparametrization invariance 
it can be fixed to an arbitrary value as long as f{t) does not equal to one. However, keeping it explicitly is quite 
convenient, since it allows to quickly derive a first-order Friedmann equation for the scale factor. 
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and the auxiliary scalars are frozen to their background values 



CT, 



X*. We define the 



perturbations of the dynamical fields in this gauge as follows, 

9,. = «'(V + K)^ ^/Mpi = In(ca) + C- (A6) 

The tensor }C, up to quadratic order in perturbations is given by the following expression 



]C^u = Si^-ca e« W — (r//^^ - h^^^ + ht^Ph.^ + ... )Sa. 



5^ - c (1 + C + ^C' + • • • )\l^^u - hP^T.^. + /i^^V^A. + . . . ) 



(A7) 



where all indices are assumed to be raised/lowered with the fiat Minkowski metric and Tj^^, = 
d'^(j)"'di,(f)'^riab = diag{(f ,1,1,1). We will also need an expansion to the quadratic order of the 
metric determinant 



_ 1 

-g = 1 + -h+ -h 

2 " 



1-2 1. . 



(A8) 



The Einstein frame action we would like to perturb can be conveniently parametrized by separating 
the pure general relativity sector in the following way 



Se = IJv^[r-qh' 



+ 



(jj HI' 



+ 3M|im2 I d^x (3 e4-/A^pi^-rfet(9^0«9^0„ . 



(A9) 



The first term describes pure GR on dS space (with a CC, consistent with the expansion rate), while 
the perturbations of the rest of the Lagrangian will describe deviation from GR. The quadratic 
perturbations of the second line of (A9) can be written as follows (note that the indices on metric 
perturbations are not raised and prime denotes a derivative w.r.t. conformal time r) 



4^) ^Ijd'xa'S^^^ {C - m?) + ^{hoo + h,.)C 



2ujH 



hoidiC + (71^00 + 72^m)C 



+ 73^00 + lihoohii + 'j5hoihoi + 'Jehijhij + 77/1 -^ + 78C^ 



(AlO) 



Here, with the definitions 

Xq = 4 + Sa-s + 4a4, xi = 2 + Sa^ + 604, X2 = 03 + 4a4, x^ = 2 + a-s + a^ , (All) 
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the coefficients appearing in the action are 



7i = ?)m^K 



[1 

[2: 
[3 
[4] 
[5] 
[6 
[7 
[8 
[9 



[10: 

[11 

[12; 
[13 

[14 

[15; 

[16 

[17; 
[is; 

[19 
[20 



72 = --m^c (9cc^X2 - 4(1 + 2c)cxi + 3(2 + c)xq) 
7s = — -n 



TTl 

74 = -^«: (A12) 

75 = TT^^ 

1 + c 

1 H^ 

7g = —m^ (3cc^X2 - 2(1 + 3c)c^xi + 3(3 + 2c)cxo - 24x3) - --(u + 6) 
16 ^4 

1 H^ 

77 = 777"^^ (2cc^xi - 3(1 + c)cxq + 12x3) + -— (w + 6) 

lo 8 

3 
7g = -m'^c (64cc^(a4 + /3) - 9c^(l + 3c)x2 + 8(1 + c)cxi - (3 + c)xo) • 
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